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The reliable detection of single quantum particles has revolutionized the field of quantum optics 
and quantum information processing. For several years, researchers have aspired to extend such 
detection possibilities to larger scale strongly correlated quantum systems, in order to record in-situ 
images of a quantum fluid in which each underlying quantum particle is detected. Here we report 
on fluorescence imaging of strongly interacting bosonic Mott insulators in an optical lattice with 
single-atom and single-site resolution. From our images, we fully reconstruct the atom distribution 
on the lattice and identify individual excitations with high fidelity. A comparison of the radial 
density and variance distributions with theory provides a precise in-situ temperature and entropy 
measurement from single images. We observe Mott-insulating plateaus with near zero entropy and 
clearly resolve the high entropy rings separating them although their width is of the order of only 
a single lattice site. Furthermore, we show how a Mott insulator melts for increasing temperatures 
due to a proliferation of local defects. Our experiments open a new avenue for the manipulation 
and analysis of strongly interacting quantum gases on a lattice, as well as for quantum information 
processing with ultracold atoms. Using the high spatial resolution, it is now possible to directly 
address individual lattice sites. One could, e.g., introduce local perturbations or access regions of 
high entropy, a crucial requirement for the implementation of novel cooling schemes for atoms on a 
lattice. 



Ultracold atoms in optical lattices have proven to be 
powerful simulators for the investigation of the static 
quantum phases and dynamical evolutions of strongly 
correlated quantum many-body systems. Prominent ex- 
amples include the quantum phase transition from a su- 
perfluid to a Mott insulator [T]-[4] and the fermionized 
Tonks-Girardeau gas for bosonic particles OE], as well 
as the recently realized fermionic Mott insulator [71 18] . In 
all these cases, the strong interactions between the parti- 
cles compared to their kinetic energy generate intriguing 
highly correlated quantum states that are of fundamental 
interest in condensed matter physics and promising for 
practical applications in quantum information science. 
For many of these applications it is of crucial impor- 
tance to image the correlated many-body systems with 
single-atom and single-site resolution. One could then 
e.g. probe the evolution from a Poissonian atom number 
distribution into highly number squeezed Fock states on a 
lattice not only globally [9 but also on a local scale. Fur- 
thermore one should be able to directly observe critical 
phenomena in the in-situ density or in spin-resolved im- 
ages of the particles when approaching a quantum critical 
point. For applications in quantum information science, 
it is essential to address and manipulate single atoms 
on individual lattice sites. A prominent example is the 
one-way quantum computer [10^ where local single par- 
ticle measurements and operations are carried out after 
a successful global entanglement operation that creates 
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FIG. 1: Experimental Setup. Two-dimensional bosonic 
quantum gases are prepared in a single 2D plane of an op- 
tical standing wave along the z— direction, which is created 
by retroreflecting a laser beam (A = 1064 nm) on the coated 
vacuum window. Additional lattice beams along the x- and 
^/-directions are used to bring the system into the strongly cor- 
related regime of a Mott insulator. The atoms are detected 
using fluorescence imaging via a high resolution microscope 
objective. Fluorescence of the atoms was induced by illumi- 
nating the quantum gas with an optical molasses that simul- 
taneously laser cools the atoms. The inset shows a section 
from a fluorescence picture of a dilute thermal cloud (points 
mark the lattice sites). 
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FIG. 2: High resolution fluorescence images of a BEC and Mott insulators. Top row: Experimentally obtained images 
of a BEC (a) and Mott insulators for increasing particle numbers (b-g) in the zero-tunneling limit. Middle row: Numerically 
reconstructed atom distribution on the lattice. The images were convoluted with the point spread function of our imaging 
system for comparison with the original images. Bottom row: Reconstructed atom number distribution. Each circle indicates 
a single atom, the points mark the lattice sites. BEC and Mis were prepared with the same in-plane harmonic confinement 
(see Supplementary Information for the Bose-Hubbard model parameters of our system). 



a highly-correlated cluster state. 

Over the past years tremendous progress has been 
made in the high-resolution and single-atom sensitive de- 
tection of atoms on a lattice [TTHTS) . However, it has 
only now become possible to apply these techniques to 
the detection of strongly correlated quantum systems, 
in the work presented here and the work of Bakr et al. 
[16 . Here we report on in-situ fluorescence imaging of a 
Mott insulator (MI) with single-atom and single-site res- 
olution. From a single image, we reconstruct the atom 
distribution on the lattice and individual thermal excita- 
tions of the MI become directly visible. This allows us 
to observe the number squeezing and the quality of an 
atomic Mott insulator down to a single lattice site. Using 
a simple model [171 US we characterize the average den- 
sity distribution and number fluctuations of the quantum 
system, and use this for an in-situ temperature measure- 
ment. We find excellent agreement between our theory 
that assumes global thermal equilibrium. Furthermore, 
we show how the incompressible Mott phase evolves into 
a compressible normal phase as the temperature is in- 
creased. 

Our experiments start with an almost pure 2D Bose- 
Einstein condensate [4 (BEC) of up to a few thousand 
^^Rb atoms that was prepared in a single pancake-shaped 
antinode of a vertical optical standing wave (beam waist 
Wo = 75 /im) oriented along the z-axis (see Fig. [T]). The 



lattice depth was Vz = 26(2) Er, where Er = / {2m\^) 
is the recoil energy, m denotes the atomic mass of ^'^Rb 
and A the lattice wavelength. Additional beams along 
the X- and ^-directions were used to load the 2D quan- 
tum gas into an optical lattice. All lattice beams had a 
wavelength A = 1064 nm resulting in a lattice period of 
532 nm. We detected the atoms in the lattice by high- 
resolution fluorescence imaging through a specially de- 
signed microscope objective with a numerical aperture 
of NA = 0.68 and an optical resolution (FWHM) of 
^ 700 nm at a wavelength of 780 nm. For detection, the 
lattice depths along all three directions were increased 
to Vx^y^z/^B ^ 300 /iK before an optical molasses in- 
duces fluorescence and simultaneously laser cooled the 
atoms l^ll^ (see Methods). In the low density thermal 
clouds (inset in Fig. [T]) individual atoms are directly vis- 
ible above an almost indiscernible background and their 
positions have a discrete spacing given by our lattice pe- 
riod (see Supplementary Information). During the imag- 
ing, atom pairs on a lattice site are immediately lost 
due to inelastic light-induced collisions [19 . We there- 
fore only detect the particle number modulo two on each 
lattice site. This essentially amounts to recording the 
parity of the atom number. 

The 2D lattice gases used in our experiments are well 
described by the Bose-Hubbard model, where particles 
are restricted to occupy the lowest energy band of the 
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lattice and their kinetic energy is characterized by a tun- 
nehng matrix element J and an on-site two-particle in- 
teraction energy U (see Refs. [20l[2T]). For a BEC loaded 
into a weak lattice potential, U/J <C 1, one expects a 
Poissonian atom number distribution on a lattice site 
as the classical coherent matter wave field of a BEC is 
characterized by Glauber's coherent states. Such states 
with an average filling of fii per lattice site, exhibit a 
corresponding variance in the particle number = n^. 
When the interactions between the particles relative to 
their kinetic energy are increased, the system under- 
goes a quantum phase transition to a Mott insulating 
state [IHS]. For homogeneous conditions and a 2D sim- 
ple square lattice, this transition is expected to occur at 
{U/J)c ^ 16.4 (see Ref. [22;), where smah shifts of this 
critical value have been reported when the system is ad- 
ditionally exposed to an underlying harmonic trapping 
potential [23j. In our case such an additional harmonic 
confinement was caused by the Gaussian beam profile of 
our lattice beams (1/e^ waist of 75/im) and resulted in an 
in-plane harmonic confinement with trapping frequencies 
a;^/(27r) = 72(4) Hz and cj^/(27r) = 83(4) Hz for lattice 
depths of 14,^ = 23(2) Er. For U/J > {U/J)c the atomic 
MI can be described by neglecting the tunneling energy 
of the system in the so called zero-tunneling approxima- 
tion. The in-trap density distribution then exhibits a pro- 
nounced shell structure of incompressible regions where 
the density is pinned to integer values and increases in a 
step-like manner from the outer wings to the inner core 
[2] |13l |24l, l25|. At zero temperature the particle num- 
ber variance at a lattice site is then expected to vanish 
(cr| = 0) resulting in perfect Fock states. For low, but 
still finite temperatures ksT <C /7, thermal fluctuations 
can be induced. These fluctuations limit the quality of 
the number squeezing and eventually lead to a complete 
melting of the characteristic shell structure of a MI when 
the temperature is increased above Tm — 0.2U/kB (see 
Refs. LlilSj. 

We monitored the dramatic differences in the density 
profiles and the on-site number fluctuations by imaging 
the in-trap atom distributions of a BEC and a MI in the 
zero-tunneling limit for different atom numbers and tem- 
peratures (see Fig.|2j first row). For the Mis, the lattices 
along the x- and ^/-directions were increased in 5-shaped 
ramps within 75 ms up to values of Vx^y = 23(2) Er. To 
freeze out the atom distribution of a BEC, we ramped up 
the lattices within 0.1ms. Using the point spread func- 
tion (PSF) of our optical imaging system we were able to 
reconstruct the atom number distribution on the lattice 
with single-site and single-atom resolution via an image 
processing algorithm (see Methods). It works reliably 
even in the regions of high atomic density, as illustrated 
in Fig.[3l 

To compare the digitally reconstructed atom distribu- 
tion (see Fig. |2j bottom row) with the raw images, we 
show the reconstructed distribution convoluted with the 
PSF in the center row of Fig. |2j For a BEC with a Pois- 
sonian atom number distribution the average filling one 




FIG. 3: Identification of single atoms in a high- 
resolution image. The points mark the centers of the lattice 
sites, circles indicate those sites where our deconvolution al- 
gorithm determined the presence of an atom. The image is a 
zoom into the upper right part of Fig. |2t. 



detects due to the parity measurement in the fluores- 
cence imaging is n^^det = 1/2 [1 — exp (— 2n^)], which sat- 
urates at n^^det = 0.5 for > 1.5. In this limit, the de- 
tected atom number variance then saturates accordingly 



at a'. 



,det 



0.25. Indeed for a BEC, we observed that 



the recorded atomic density exhibits large atom number 
fluctuations from site to site. In contrast, for a MI we 
obtain plateaus of constant integer density, with almost 
vanishing fluctuations. For increasing particle numbers, 
the images in Fig. [2] show how successive MI shells are 
formed, which appear as alternating rings of one and zero 
atoms per site due to our parity measurement. In both 
the raw images and the reconstructed ones, individual 
defects are directly visible. The high symmetry of our 
atom clouds reflects the high optical quality of our lattice 
potentials. A small ellipticity is caused by the different 
harmonic trapping frequencies ujx and ujy. 

We used the reconstructed site occupation numbers to 
determine the temperature of the sample based on a sin- 
gle image. For deep lattices, U / J 300, as used in our 
experiments for Mis, tunneling becomes completely sup- 
pressed such that coherent particle-hole fluctuations are 
expected to be negligible and defects are only induced 
by thermal fluctuations. The symmetry of our clouds 
allowed us to average the data azimuthally, taking into 
account the ellipticity, and to obtain radial profiles for the 
average density ndet(^) and variance crdet(^) (^^^ Fig.|4^,b 
and Methods). We fitted analytic expressions derived in 
the zero-tunneling regime (see Methods) to our data. The 
results of such a fit for an n = 1 (0 < fi/U < 1) and an 
n = 2 (1 < iJi/U < 2) MI are displayed in Fig. |4^,b. 
The MI regions can be identified as connected regions 
of constant integer density and vanishing on-site number 
fluctuations, which in the atomic limit of the Hubbard 
model signify the presence of incompressible Mott do- 
mains [23 . For both density profiles and atom number 
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FIG. 4: Radial atom density and variance profiles. Ra- 
dial profiles were obtained from the digitized reconstructed 
images by azimuthal averaging, a, b, Yellow and red points 
correspond to the n = 1 and n = 2 MI images of Fig.[2]l,e. 
The grey points were obtained from a BEG (data from 
Fig. [2^) for reference. The displayed statistical errors are 
Glopper-Pearson 68% confidence intervals for the binomially 
distributed number of excitations. For the Mis both den- 
sity and variance profiles are fitted simultaneously with the 
model functions of Eqs. ([T]) and ([2| (see Methods) with T, 
fi and ro = 2U / (jnuJxOOy) as free parameters. For the two 
curves, the fits yielded temperatures T — 0.090(5)[//A:s and 
T = 0.074(5)[//A:s, chemical potentials /i = 0.73(3)[/ and 
jjL — 1.17(1)[/, and radii ro = 5.7(1) /xm and ro = 5.95(4) /xm 
respectively. From the fitted values of T, /i and ro, we de- 
termined the atom numbers of the system to N — 300(20) 
and N — 610(20). c, d. The same data plotted versus the lo- 
cal chemical potential using the local-density approximation. 
The inset of c is a Bose-Hubbard phase diagram (T = 0) 
showing the transition between the characteristic MI lobes 
and the superfluid region. The line starting at the maximum 
chemical potential /i shows the part of the phase diagram ex- 
isting simultaneously at different radii in the trap due to the 
external harmonic confinement. The inset of d is the entropy 
density calculated for the displayed n = 2 MI. 



variances we find excellent agreement between the ex- 
perimental data and the theoretical model for all radial 
distances. This supports the assumption that our system 



is in global thermal equilibrium, in contrast to Ref. [26] . 
The extracted temperatures of T = 0.090(5)[///cb and 
T = 0.074(5) for the n = 1 and n = 2 data are 
well below the MI melting temperature T^. Our tem- 
perature estimates are conservative since all defects are 
attributed to thermal excitations in our model. How- 
ever, defects might also arise due to "collateral damage" 
caused by atoms undergoing the light-induced collisions. 
For reference, we show the corresponding data obtained 
by freezing out the atom distribution of a BEG. We ob- 
serve the expected saturation of ndet at 0.5 together with 
a maximum variance of a^^^ at 0.25. We note that the 
thermal shells surrounding a MI core also exhibit this 
maximum variance and can be as narrow as 1-2 lattice 
sites. In Fig. we plot both MI data sets versus local 
chemical potential. In a single image, we thus mapped 
out an entire line in the phase diagram as illustrated in 
the inset of Fig. |4j3. The slightly different temperatures 
of the two Mis are clearly visible in the different widths 
of the variance curves. 

Our measurements also confirm with unprecedented 
clarity that the entropy of the strongly correlated quan- 
tum gas is concentrated around the Mott insulating re- 
gions, whereas in the center of a MI, for local chemi- 
cal potentials of /iioc = (n + l/2)/7, number fluctuations 
are completely suppressed and the entropy density is es- 
sentially zero. For the lowest observed temperature of 
T = 0.074{b)U/kB we calculate a 99.7(1)% probability 
for unity occupation in the center of the n = 1 Mott- 
insulating plateau. Using the zero-tunneling model we 
can also extract the total entropy per particle for our sys- 
tem S/{NkB) = 0.34(2) (see Methods) which is around 
the critical entropy for quantum magnetism [27 . 

Finally, we show how a Mott insulator melts, as the 
temperature (or entropy) of the quantum gas is increased 
(see Fig.[5|. At constant total chemical potential but 
increasing temperatures, one observes that the Mott do- 
mains gradually vanish. Although there is no sharp tran- 
sition to a normal fluid state in this case, Mott plateaus 
and number squeezing degrade rapidly, once T > Tm (see 
Figs.[5]i,e), as predicted in Ref. [17]. 

In summary, we have demonstrated single-site and 
single-particle resolved detection of a strongly corre- 
lated system in an optical lattice. Our method can 
be extended to investigate quantum critical phenom- 
ena, density-density correlations or even non-local string- 
operators that are inaccessible in condensed matter ex- 
periments. Our imaging system can be used to focus an 
off-resonant laser beam onto a single atom and thereby 
allow single-site manipulation of the atomic qubits [28^. 
This opens a new avenue for experiments with ultracold 
quantum gases, where e.g. novel cooling schemes may 
be applied by accessing regions of high entropy [29] . For 
future work it would be interesting to investigate how 
entropy propagates in strongly correlated systems, after 
injecting it on a local scale into the system. Atoms in the 
MI with one atom per lattice site are also very promising 
as a quantum register for scalable quantum computing. 
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FIG. 5: Melting of a Mott Insulator, a-c Strongly 
correlated atomic samples for three different temperatures, 
and constant total chemical potentials in the zero-tunneling 
limit. For higher temperatures, an increased number of inde- 
pendent particles or holes appear. Data shown was binned 
2x2. e, f. Density and variance profiles as a function of 
chemical potential, determined as described in the caption 
of Fig.|4] Red, orange and yellow points correspond to the 
datasets from a, b, and c, respectively. Light blue points 
correspond to the low-temperature n = 2 MI of Fig. [2^ and 
Fig.|4]with T = 0.074{5)U/kB and /i = 1.17(1)[/. The pa- 
rameters extracted from the radial fits are T = 0.17(1)[//A;s, 
/i = 2.08(4)[//A:s, for a, T = 0.20(2)[//A:s, M = 2.10(5)[//A:s, 
for b and T = 0.25(2)[///cs, M = 2.06(7)[///cs, for c. 



especially with the very low defect density shown in this 
work. 
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Methods 

Preparation of a single 2D system. We started 
by loading atoms from a magnetic quadrupole trap into 
a single beam optical dipole trap (1064 nm, beam waist 
wq = 40 /im). By translating the focus of the dipole 
trap using mirrors on a motorized micrometer stage, the 
atoms were transported in front of the high resolution 
imaging system within 2.5 s. A magnetic quadrupole field 
whose center is shifted below the trap compressed the 
cloud in the axial direction of the dipole trap laser beam. 
After 500 ms of evaporative cooling in this hybrid trap 
configuration, we transferred the atoms into the z-lattice. 
It is oriented along the optical axis of the imaging sys- 



tem, parallel to gravity, and superimposed with the hy- 
brid trap. Initially we populated 60 antinodes (slices) of 
the standing wave, creating independent 2D systems. To 
extract a single 2D system, we used position dependent 
microwave transfer in a magnetic field gradient. The gra- 
dient with dB/dz = 24G/cm (together with a bias field 
of 32 G) was produced with a single coil placed 50 mm 
above the atoms; the coil axis coincides with the z-lattice 
beam. Our magnetic field gradient results in a position 
dependent frequency shift dv/dz ^ 5 kHz//im of the 
\F = l^rriF = —1) ^ \F = 2,mi? = —2) transition. An 
initial microwave frequency sweep over 1 MHz brought 
all atoms from \F = I^tyif = —1) to \F = 2,mi? = —2). 
We then transferred atoms from one slice back to the 
|F = 1, rriF = — 1) state using a resonant Blackman pulse 
of 5 ms duration. All atoms remaining in F = 2 were re- 
moved from the trap by a laser pulse resonant with the 
F = 2 ^ F' = 3 transition. We then evaporatively 
cooled the atoms by ramping down the intensity of the 
z-lattice from 300 Er to ~ 25 Er within 1.5 s, while si- 
multaneously tilting the potential along the horizontal 
direction with a magnetic field gradient [30 . Depending 
on the end point of this evaporation, we created BECs 
with atom numbers ranging from 50 — 2000. 

Imaging single atoms in the lattice. Our micro- 
scope objective was custom made (Leica Microsystems) 
and is located outside of the vacuum chamber with a 
working distance of 13 mm. We detected the atoms by 
illuminating them with an optical molasses, red detuned 
with respect to the free space resonance by 45 MHz. It 
consists of two pairs of retroreflected laser beams su- 
perimposed with the horizontal lattice axes. A third z- 
molasses laser beam aligned in reverse direction through 
the imaging system provided cooling in the vertical di- 
rection. The total scattering rate from all laser beams 
was ~ 60 kHz. With our total detection efficiency of 

9 % (solid angle 15%, transmission of all optical el- 
ements 71%, camera quantum efficiency 85%), we col- 
lect about 5000 photons/atom within our illumination 
time of 900 ms. To ensure a homogeneous illumination 
of all atoms, we scanned the retro-reflecting mirrors of 
the horizontal molasses beams (at a frequency of 300 Hz 
and 400 Hz) using piezo elements and mutually detuned 
the molasses beams by 43 Hz. Additionally, we spatially 
scanned the z-molasses beam across the cloud at a fre- 
quency of 100 Hz and an amplitude of 35 /im. We also 
corrected for an etaloning effect of the CCD camera, 
which caused a spatially dependent signal strength. We 
carefully optimized the molasses parameters to minimize 
hopping of the atoms to adjacent lattice sites, by taking 
two consecutive images of the same cloud for the first 
400 ms and the last 400 ms of our 900 ms illumination 
period. From the analysis of several of these double im- 
ages, we found that about 0.5% of the atoms hop during 
the molasses illumination. Before switching on the mo- 
lasses, we removed the atoms in the doubly occupied sites 
with a 50 ms pulse on the F = 2 to F' = 3 transition, 
which is 6.8 GHz red detuned for the atoms in F = 1, 
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but efficiently excites into the molecular potentials. 

For the horizontal optical lattice laser beams we used 
two single mode fiber amplifiers seeded with the same 
narrowband solid-state laser, whereas the vertical lattice 
beam was derived from an independent solid-state laser. 
The horizontal axes were mutually detuned by 220 MHz 
and had orthogonal polarizations. We obtained about 
low per lattice axis at the experiment, yielding trap 
depths of about 300 fiK. 

Reconstruction of the atom number distribu- 
tion. We developed a deconvolution algorithm to re- 
construct the atom number distribution from a fiuores- 
cence image. It uses a model of the point spread func- 
tion (PSF) of our imaging system that was determined 
from averaging over many images of isolated individual 
atoms (see Supplementary Information). The algorithm 
tries different model-configurations for each lattice site 
and its nearest neighbors in order to minimize the dif- 
ference of the original image with the reconstructed one. 
This reconstructed image was obtained by convoluting 
the atom number distribution with our PSF (see centre 
row of Fig. [2|. The algorithm allows for a variance of 
the fiuorescence level of each atom within ±20% of the 
mean photon counts. These varying fiuorescence levels 
partially arise from the inhomogeneous intensity of the 
molasses light. We additionally found an increased fiu- 
orescence level of about 10%-20% in the center of very 
dense n = 1 shells of a MI, compared with the isolated 
atoms in the outer part of the images. This effect might 
arise from the partial coherence of the light scattered by 
the atoms, combined with their regular distribution in 
the lattice. 

We have evaluated the fidelity of the reconstruction al- 
gorithm by creating simulated images of a known atom 
distribution using the PSF of our imaging system, the 
poissonian and superpoissonian noise contributions of the 
light hitting the EMCCD camera (including the ampli- 
fication process), and the site-to-site fiuorescence fiuctu- 
ations of ±20%. Running the reconstruction algorithm 
over several hundred of such randomly generated images 
of Mis at finite temperatures, we find a reconstruction 
fidelity of ~ 99.5%. In our experiment, the main limita- 
tions of the fidelity are atom losses during the detection 
process due to collisions with background gas atoms. We 
measured that about 1% of the atoms are lost during 



the 900 ms detection period, which corresponds to a trap 
lifetime of ~ 75 s. 

Radial atom number distribution and variance. 

In the zero-tunneling regime (the "atomic limit" of a 
MI), the atom number distribution at a lattice site at 
radius r is given by P^(n) = e^[^i°'=^^^^~^"']/Z(r), where 
^(^) — e^[^^°'=^^^^~^'^] is the grand canonical par- 
tition function, /3 = l/(/c^T), //loc(^) the local chemi- 
cal potential and En = Un{n — l)/2 is the interaction 
energy. Using a local density approximation (LDA), 
we define /iioc in terms of the global chemical poten- 
tial /i and the external harmonic trapping confinement: 
Mioc(^) = ^^(^x^^ ±^y^^)- Taking the light induced 
losses into account we calculate the expected detected 
density at different radii: 

ndet(r) = ^ E mod2(n)e''[''-('-)"-^"l (1) 

In the presence of light induced collisions ^^et(^) ~ 
^det(^) and the detected variance is therefore simply 

^det (^) = ^det (r) - nl^t (r) ' (2) 

We extracted radial density and variance profiles from 
the reconstructed two-dimensional atom distribution of 
a single image. For this, we first determined the center 
of the cloud, and then binned the lattice sites according 
to their distance from the center, thereby correcting for 
the ellipticity of 10%. The bin sizes were chosen larger 
near the center to have sufiicient statistics. 

We fitted the experimental profiles to Eqns. ([T]) and 
([2| and extracted the temperature and the global chem- 
ical potential. These can then be used to calculate 
the original atom number distribution P^(n). Insert- 
ing the radius corresponding to /iioc = 0.5 we extract 
the maximal theoretical unity occupation probability. 
We can furthermore calculate the local entropy density 
5'ioc(^) = —kB'^n-^r{'n)l^[Pr{'n)]. Summing the density 
and entropy density over the lattice sites we calculate the 
total number of particles N = 591(9) and the total en- 
tropy S/kB = 200(13) given the fitted values of T and fi 
from the n = 2 data of Fig. |4] This gives the entropy 
per particle S/iNks) = 0.34(2). 
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Supplementary Information 

Determination of lattice angles and spacing. To 

characterize our imaging system and to determine the 
lattice structure, we used a fluorescence image of a dilute 
thermal could, similar to the inset in Fig. [l] of the main 
text. The lattice axes are oriented at approximately ±45° 
with respect to our images. A precise determination of 
this angle and the lattice spacing is needed so that the 
deconvolution algorithm works with high fidelity We 
first determined the center positions of isolated atoms 
from this image by a simple fitting algorithm. The his- 
togram of the mutual distances projected in a coordinate 
system rotated by an angle clearly shows the periodic- 
ity of the lattice (see Fig. [6|i,b) and the visibility of the 
pattern depends very sensitively on 6. For a quantitative 
analysis, we fit a sum of equidistant Gaussians to the his- 
togram. The width of the Gaussians for different values 
of (Fig. ^) shows a clear minimum at = 45.85(1)°. 
We obtained a similar graph for the other lattice axis and 
found an angle of —45.55(1)°. The distance of the Gaus- 
sians is 4.269(4) pixel which corresponds to lattice period 
of 532 nm. Thus, our magnification factor is 128.4(2) and 




Coordinate rotation angle (°) 



FIG. 6: Determination of the lattice angles. a,b His- 
togram of the distances between the center positions of indi- 
vidual atoms projected in a coordinate system rotated by an 
angle 0. The line is a fit to a sum of equidistant Gaussians. 
c The width of the fitted Gaussians show a clear minimum. 
The red line is a parabolic fit and yields a minimum rotation 
angle at 6* = 45.85(1)°. 
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FIG. 7: Azimuthal average of our experimentally ob- 
tained point spread function. The data was obtained by 
averaging over 68 signals of single atoms. The line is a fit 
with the double Gaussian of Eq. (|3} and yields ai = 2.06(5) 
pixels, (72 = 9.6(1.2) pixels, a = 0.075(2) and C = 1050(7). 

one pixel of the CCD camera corresponds to 124.6(1) nm 
in the object plane. The angles and lattice spacing de- 
termined by this method are used as fixed parameters 
for our deconvolution algorithm. We also found that the 
phases of the two lattice axes slightly drift from shot to 
shot. They are determined for each image by fitting the 
center positions of single atoms in the outer part of the 
images. 

Determination of the point spread function. We 

determined our point spread function (PSF) from the flu- 
orescence image of a dilute atomic cloud. We summed 
the fluorescence image of many individual atoms that 
were isolated from their neighbors by more that 12 pix- 



els. The summed image is almost radially symmetric and 
we computed an azimuthal average (see Fig. [t]). We ex- 
pect our PSF to be a convolution of an Airy disk with a 
Gaussian, taking into account fluctuations of the lattice 
with respect to the imaging system and the width of the 
atomic wavepacket in the potential wells. Due to this 
convolution, the first minimum of the airy pattern is not 
visible in our averaged signal. We found that our PSF 
can be well approximated by a double Gaussian: 

PSF{x,y) = C[{l-a)exp{-0.5{x^ ^y^)/al) 

+ aexp(-0.5(x2+^2)/ai)] (3) 

with widths ai, (72 and a parameter a describing the 
relative amplitudes. The maximum fluorescence level C 
varies from day to day and is in the range of 800-1200 
counts. 

Single-band Hubbard parameters 

For a z-lattice depth of Vz = 26Er^ the single-band 
Hubbard parameters for our 2D system are given in the 
following table: 



Vx,y{Er) 


J/h (Hz) 


U/h (Hz) 


5 


134 


382 


10 


39 


607 


15 


13 


776 


20 


5 


917 


25 


2 


1039 



Note that one expects the the single-band Hubbard pa- 
rameters to be slightly renormalized due to multi-orbital 
effects. 



